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Abstract. We study, within an effective field theory framework, 0(E 2 /Mpy) Planck- 
i' scale suppressed Lorentz invariance violation (LV) effects in the neutrino sector, whose 

size we parameterize by a dimensionless parameter r/ u . We find deviations from 
predictions of Lorentz invariant physics in the cosmogenic neutrino spectrum. For 
positive 0(1) coefficients no neutrino will survive above 10 19 eV. The existence 
of this cutoff generates a bump in the neutrino spectrum at energies of 10 17 eV. 
Although at present no constraint can be cast, as current experiments do not have 
OA ■ enough sensitivity to detect ultra-high-energy neutrinos, we show that experiments in 

construction or being planned have the potential to cast limits as strong as r\ v < 10~ 4 
on the neutrino LV parameter, depending on how LV is distributed among neutrino 
mass states. Constraints on r\ v < can in principle be obtained with this strategy, but 
Q\ ' they require a more detailed modeling of how LV affects the neutrino sector. 

o ' 



1. Introduction 

Over the last fifteen years there has been consistent theoretical interest in possible small 
deviations from the exact local Lorentz Invariance (LI) of general relativity as well as a 
flourishing of observational tests. The theoretical interest is driven primarily by ideas 
in the Quantum Gravity (QG) community that Lorentz invariance may not be an exact 
local symmetry of the vacuum. The possibility of outright Lorentz symmetry violation 
(LV) or a different realization of the symmetry than in special relativity has arisen in 
string theory [TJ |2], Loop QG [31 IH E], non- commutative geometry [El El El E], space- 
time foam [10], some brane-world backgrounds [TT] . and condensed matter analogues 
of "emergent gravity" [12] . As well, allowing the fundamental theory of gravity to be 
non-relativistic in the ultraviolet can make gravity renormalizable while avoiding some 
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of the other pathologies that plague renormalizable gravitational actions with higher 
derivative terms [T3] . 

Constructing useful tests for the various active models and ideas is therefore vital. 
Since there are so many theoretical models around, a good approach is to work within 
a calculable framework where all possible LV terms are parameterized. Each theory 
then picks out a certain combination of terms which can be constrained. In this vein, 
a standard method is to simply analyze a Lagrangian containing the standard model 
fields and all LV operators of interest that can be constructed by coupling the standard 
model fields to new LV tensor fields that have non-zero vacuum expectation valued- 

All renormalizable LV operators that can be added to the standard model in this 
way are known as the Standard Model Extension (SME) [TB]. These operators all have 
mass dimension three or four and can be further classified by their behavior under CPT. 
Higher mass dimension operators can be systematically explored as well, which is useful 
in case the naive EFT hierarchy breaks down due to other new physics (e.g. SUSY) or 
quantum gravity introducing a custodial mechanism for the renormalizable operators in 
the infrared. 

The CPT odd dimension five kinetic terms for QED coupled to a non-zero vector 
field were written down in [17] while the full set of dimension five operators with a 
vector were analyzed in [TS]. The dimension five and six CPT even kinetic terms for 
QED for particles coupled to a non-zero background vector, which we are primarily 
interested in here, were partially analyzed in [19] ■ The full set of dimension five and six 
operators for QED has recently been introduced in [20] . It is notable that SUSY forbids 
renormalizable operators for matter coupled to non-zero vectors [21] but permits certain 
nonrenormalizable operators at mass dimension five and six. 

Many of the parameterized LV operators have been very tightly constrained via 
direct observations (see [221 [231 EEH [21] for reviews). In particular, the dimension 
five and six CPT even operators for LV with a vector field have been recently 
directl}j§| constrained in the hadronic sector by exploiting ultra-high energy cosmic- 
ray observations [26] performed by the Pierre Auger Observatory (PAO). Indeed, the 
construction and successful operation of this instrument has brought UHECRs to the 
interest of a wide community of scientists and it is expected to allow, in the near future, 
the assessment of several problems of UHECR physics and also to test fundamental 
physics (in particular Lorentz invariance in the QED sector) with unprecedented 

| There are other approaches to either violate or modify Lorentz invariance, that do not necessarily 
yield a low energy EFT (see [T3] and refs therein). However, these models do not easily lend themselves 
to particle physics constraints as the dynamics of particles is less well understood and hence we do not 
consider them here. In particular, we remark here that ideas of deformation, rather than breaking, of 
the Lorentz symmetry (see, e.g., [TS]) do not have an ordinary-EFT formulation, hence they cannot be 
tested with the arguments presented in this work. 

§ Notice that all these operators can be indirectly constrained by EFT arguments [55], as higher 
dimension LV operators induce large renormalizable ones if we assume no other relevant physics enters 
between the TeV and Mp\ energies. SUSY, however, is an example of new relevant physics that can 
change this. 
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precision [23 EEl EE] ■ 

The UHECR constraints [301 EH E21 E31 EH ESI ESI EH] rely on the behavior 
of particle reaction thresholds with LV, which are one of the best methods in the 
EFT approach to constrain nonrenormalizable LV operators. Many LV operators 
give modified dispersion relations for free particles, where the energy as a function of 
momentum deviates slightly from the special relativistic form. For threshold reactions 
what matters is not the size of the LV correction to the energy compared to the absolute 
energy of the particle, but instead the size of the LV correction to the mass of the 
particles in the reaction. Hence the LV terms usually become important when their size 
becomes comparable to the mass of the heaviest particle. If the LV term scales with 
energy as E n , then this critical energy is E cr ~ (m 2 M"f 2 ) 1/n [37J. According to this 
reasoning, the larger the particle mass the higher is the energy at which threshold LV 
effects come into play. This is why > TeV electrons and positrons, but not protons, 
can be used to constrain n = 3 LV [3S], and why UHE protons are needed to obtain 
constraints on hadronic LV with n = 4 scaling (which corresponds to CPT even mass 
dimension five and six operators). 

From this point of view, neutrinos, with their tiny mass of order m v ~ 0.01 eV 
|39j, are in principle the most suited particles to provide strong constraints on LV, 
at least for reactions involving only neutrinos. One such reaction is that of neutrino 
oscillation. Indeed, for a decade neutrino oscillations have proven to be excellent tests 
of the SME and other LV models [HI1 EH HD E21 S3] as when the LV corrections are 
near the neutrino mass, the oscillation pattern can change as a function of energy, 
direction and mass. ICECUBE may even be able to probe dimension six operators 
with time of flight techniques with TeV neutrinos from distant Gamma Ray Bursts [44]. 
UHECR experiments have also the capability of placing constraints on SME parameters 
by exploiting neutrino oscillations [45|jjJ. One can construct complementary and in some 
measures even more sensitive neutrino tests of higher dimension operators by leveraging 
observables which deviate more strongly from their special relativistic values as the 
neutrino travel distance increases. One such observable is found to be related to UHE 
neutrino spectrum observations. 

Despite the threshold being low for LV effects to kick in, neutrinos with ultra-high 
energy are necessary to achieve a signal, as they interact so weakly that the phase 
space for a LV reaction must be huge to generate an appreciable rate. This requirement 
implies that, since the LV terms and hence the phase space grow with energy, very 
large energies are needed. Indeed, for renormalizable operators, where the phase space 
does not grow quickly enough, reactions of the type we consider here never achieve 
the necessary rate given current bounds on their coefficients, even for UHE neutrinos 
over cosmological distances [451 I4D] . However, for the CPT even non-renormalizable 
operators the situation is different, and we find that for energies of order 10 17 eV there 
can be a significant modification of the spectrum. Next generation neutrino detectors 

| We notice that in the same work [45] the process of neutrino decay is discussed, but in a different 
context than what we consider here. 
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such as ANITA [UJ and SuperEUSO [HI H9] are sensitive to neutrinos of energies 
> 10 19 eV. Further experiments, like the planned ARIANNA [SUl EE] and IceRay [52J, 
will cover the range 10 17 -r- 10 20 eV. 

In the present work we study how limits on the absolute scale of non-renormalizable 
CPT even LV neutrino parameters can be obtained from UHE neutrino observations. In 
particular we will determine the energy scale where the neutrino spectrum might begin 
to deviate as a function of the size of the LV coefficients. With three neutrino species 
as well as assorted light leptons possibly involved, the spectrum as a function of various 
LV parameters can only be computed by a detailed parameterized numerical search. 
At this stage, where LV in this sector is only speculative, we feel that such a search is 
unwarranted. However, we will show the "best case" scenario for a LV neutrino signal 
with dimension six operators as well as a scenario where a signal in the UHE spectrum 
is much more subtle even though LV is still relatively strong. 

This paper is structured as follows. In section [2] we describe the theoretical LV 
framework in which we will derive the LV effects on the UHE neutrino spectrum. In 
sections [3] and H] we will give general information on the standard understanding of 
UHE neutrino generation by UHECRs, and we will present and detail the main LV 
reactions possibly affecting their spectrum. Section [5] is devoted to present results in 
our test cases. In Section Owe discuss the possible role of other processes than neutrino 
splitting. In section [7J we will report our final remarks and conclusions. 



2. Theoretical framework 



In order to study the phenomenological consequences of LV induced by QG, the existence 
of a dynamical framework in which to compute reaction thresholds and rates is essential. 
We assume that the low energy effects of LV can be parameterized in terms of a local 
EFT. Furthermore, for simplicity we assume that only boost invariance is broken, while 
rotations are preserved (see [22] for further comments on rotation breaking in this 
context). Therefore we introduce LV by coupling standard model fields to a non-zero 
vector. More complicated forms could of course be chosen, however this would introduce 
significant direction dependence. With low statistics already for neutrino observatories, 
disentangling direction dependence would be a difficult task. 

We focus on the CPT even mass dimension five and six operators involving a vector 
field u a (which we assume to describe the preferred reference frame in which the CMB is 
seen as isotropic) coupled to a Dirac neutrino ip in a mass eigenstate with mass m that 
are quadratic in matter fields and hence modify the free field equations. Neglecting the 
left-handed helicity of neutrinos for a moment, the Lagrangian for a generic fermion is 
the usual Dirac term plus [19| 

1> 



-!-{ u .Df{afP L + afP R ) (1) 
{u-Dfiu-^iafPL + afPa) 
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M P1 



E?=p t + m* + f£>p l + f£>jj T , (2) 



where u a is a timelike unit vector describing the preferred frame, Pr and Pl are the 
usual right and left projection operators, Pr,l = (l±7 5 )/2, and D is the gauge covariant 
derivative. The a coefficients are dimensionless. 

For fermions, at E 3> m the helicity eigenstates are almost chiral, with mixing due 
to the particle mass and the dimension five operators. Since we will be interested in 
high energy states, we re-label the a coefficients by helicity, i.e. = , = otf* . 
The resulting high energy dispersion relation for positive and negative helicity particles 
can easily be seen from (TjQ) to involve only the appropriate at+ or a_ terms. For 
compactness, we denote the helicity based dispersion by a± . Therefore at high energies 
we have the dispersion relation 

p2 _ J2, L _2 _l ^(4)^2 , V 

Z P1 

where = j^{ol^ + a+ ) and f± ^ = 2a± + a^a+K We have dropped the 
terms as the □ operator present in these terms makes the correction to the equations 
of motion proportional to m 2 and hence tiny. 

The dimension five fermion operators induce two corrections, one proportional to 
E 4 and one corresponding to a change in the limiting speed of the fermion away from 
c. The second effect, generated by is naturally of order 10~ 30 due to the mass 
suppression for neutrinos and hence can also be disregarded. Furthermore, we can 
drop the positive helicity coefficients, as the dominant signal will be from left-handed 
neutrinos produced by standard model couplings. Finally, since CPT is conserved, and 
neutrinos and antineutrinos exist only in opposite states of helicity, = r\ v = fjf\ 
All these lovely properties of neutrinos mean we have just one LV coefficient, while in 
other cases [26J further assumptions, such as parity preservation, are needed to reduce 
the number of free parameters. The final dispersion relation we assume in this work for 
each mass eigenstate neutrinos and antineutrinos is 

4 

E 2 U =p 2 + ml + r] uI — T , (3) 

where / denotes the mass eigenstate of the neutrino. There is no reason why each mass 
eigenstate needs have the same coefficient, and indeed whether or not they do makes a 
dramatic difference in the observed spectrum. In particular, we will study two cases: in 
the "flavor blind" case all the r\ v i have roughly the same magnitude, which translates 
in all the LV effects being essentially independent of flavor; in the "flavor dependent" 
case the r] u j are instead different for different mass states, which makes the LV effects 
depend upon flavor states. Among the "flavor blind" cases we can find a "best case" 
scenario, in which the effects of LV are maximal, as well as a "worst case" scenario, in 
which LV is present but is ineffective in the context of UHE neutrinos we are discussing 
now. 

A fairly accurate general estimate of the minimum energy in which LV corrections in 
equation (j^J) is relevant is obtained, as stated previously, by comparing the largest mass 
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of the particles entering in the LV reaction with the magnitude of the LV correction 
in these equations [37J. In our case, assuming t] v i ~ 1, the typical energy at which 
LV contributions start to be relevant is of order E t h ~ \/m u Mp\ ~ 10 13 eV, assuming 
m v ~ 10~ 2 eV. This energy of 10 TeV allows observatories such as ICECUBE to possibly 
constrain these operators in the future using neutrino oscillations. We instead focus on 
changes to the UHE neutrino spectrum induced by these operators. 

3. UHE neutrinos and LV 

If we neglect exotic sources of UHE neutrinos (as suggested in many top-down models 
for the production of UHECRs, now disfavored by the current experimental photon 
limits), the "cosmogenic" neutrino flux is created [531 E2 E3 EH] via the decay of 
charged pions produced by the interaction of primary nucleons with CMB photons above 
E w - 5x 10 19 eV, the Greisen-Zatsepin-Kuz'min (GZK) effect [57]. HiRes [5g] and 
AUGER [59] spectral observations seem to confirm the presence of a GZK suppression 
in the UHECR spectrum. Although the suppression of the UHECR spectrum could be 
also due to the maximal accelerating power of UHECR sources, the fact that it occurs 
at just the right energy for being GZK taking place during propagation, and the results 
|60] on the correlation of the UHECR arrival directions with the large scale distribution 
of matter within ~ 75 Mpc, seem to favor the GZK explanation. Motivated by these 
considerations, we make here the hypothesis that the GZK reaction is at work during 
UHECR propagation. 

We have checked numerically that the process of pion decay is not strongly affected 
by the above LV if only neutrinos are LV. One might ask in general if LV effects in 
the hadronic sector, which we have so far neglected, can matter in the production of 
UHE neutrinos. We can neglect the hadronic sector because of the tight constraints 
already placed on such LV operators [271 E3 E3 EE] , which are stronger than what will 
be considered here. 

Violation of Lorentz invariance however introduces new phenomena in the 
propagation of UHE neutrinos. A partial list of these effects includes: 

Modified v oscillations Since the LV is parameterized in the mass eigenstate, the 
LV terms act as contributions to an effective mass and contribute to neutrino 
oscillations. Cosmogenic neutrinos are not the right phenomena with which to 
study modified oscillations as they have oscillated many times during their flight 
from an unknown source, making it extremely hard to derive oscillation constraints. 
ICECUBE and other neutrino observatories may be sensitive to these operators for 
atmospheric neutrinos. 

z/-Cerenkov Cerenkov radiation in vacuum via charge-radius coupling or gravity: 
v — > uj. This possibility has been already investigated for renormalizable 
operators [1H1 HQ]. The rate is too small for both renormalizable operators (at 
current limits) or non-renormalizable operators, even for cosmogenic neutrinos. It 
also involves additional LV coefficients for photon emission and perhaps new modes 
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for graviton emission [61], thereby complicating the analysis and, at best, yielding 
less definitive limits. 

//-splitting : v — > vvv. This effect is what we focus on in the following, as it exclusively 
involves the neutrino sector and has a high enough rate to be seen at UH energies. 

z/-pair emission : v — > vZ — > vff. Here, / represents some other fermion species 
besides neutrinos. Electrons are the only fermion species light enough to have an 
effect compared to //-splitting at UH energies. We initially will ignore this reaction, 
but we will consider the possible effects of electron pair and hadronic emission in 
section [HI 



4. Neutrino Splitting 

To understand how a neutrino can split into three, it suffices to calculate the threshold 
energy for this reaction to occur, which in a Lorentz invariant scenario would be infinite. 
Let us compute the threshold for the decay ua{p) j/ a(p') 1 'b(q)^b(q') , where z/4, v B are 
neutrinos of mass m^,ms and p, p', q, q' are the momenta. In the threshold configuration 
the momenta of the outgoing particles are all aligned and parallel to the direction of the 
incoming neutrino momentum [62J, hence we choose all momenta to be in the z-direction. 
Let x,y and t be the fraction of initial momentum carried respectively by the outgoing 
z/4, by Vbi and by z/g, t = 1 — x — y and < x, y, t < 1. For a general dispersion relation 
with LV term r]^p n /Mp^ 2 the threshold equation (the energy conservation equation 
imposing momentum conservation) can then be written as 

jij, foffa - x"- 1 ) - *.(»"- 1 + 1- 1 )] = <^+< (I + ]) -(4) 

When LI is exact, the left hand side vanishes and there is no solution. With LV, in 
the simple case where n = 4 and A = B, we obtain the equation 



P 4 (4) = _}__ 

M v\ m l A %A 3x yt 



€ =7T-. (5) 



L pr' b u A 

The minimum of the right-hand side in the equation above is 9 and is attained at 
x = y = t = 1/3, hence the threshold energy is given by 



v—rvv 
Pth(4) 



3Mpim ;; 20 TeV 



mvA (6) 



(4) 



1/2 



Note that r\ VA must be positive in order for there to be splitting, i.e. we must 
deal with superluminal neutrinos. If A 7^ B then the threshold energy must be solved 
for numerically as the LV parameter and mass are generically different. However, the 
existence of a finite threshold can still be shown for appropriate values of rj UA and rj UB . In 
the latter case, the coefficients do not need to be all positive for there to be a threshold, 
however it is necessary that r\ VB < r\ VA — K(p), where K{p) is a positive number that 
depends on the incoming momentum. Thresholds for negative coefficients have been 
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investigated previously [37J E3] and there are strong limits on the LV coefficients^ 
However, for UHE neutrinos it is more difficult to extract constraints, as we shall discuss 
in Section 15721 



4-1. The best case scenario: flavor blind LV 

Neutrinos created in a flavor eigenstate are, of course, a combination of mass eigenstates. 
In order to maximize the neutrino splitting signal one would want all mass eigenstates 
to split ("flavor blind" scenario). The best case scenario is therefore when r\ v i > for 
each / and each rj uI is of the same order of magnitude (so all states begin to decay at 
the same energy). As our best case scenario then we take r\ v i equal and positive for all 
I. 



4-1.1. Decay time computation Let us now consider, within the above "flavor blind 
scenario", energies well above threshold so that neutrinos do effectively split, and the 
best case combination of the r) u i- The splitting rate can be calculated from the neutrino 
width 

Ftot,A = F AB , (7) 
B 

where the sum runs over the open splitting channels and Tab represents the partial 
width for the channel z/4 — > vaVbVb- Tab is simply 

1 



AB 



2E P 



x / \Mab\ 2 x d^ AB , (8) 



where $ab represents the phase space of the final states and M.ab the matrix element 
for the process. The dominant channel for neutrino splitting is via the tree level neutral 
current interaction. The matrix element for this interaction is 

Mab oc —^-u A {p')^P L u A {p) - JT u B (q')rP L VB(q) , (9) 

4 cos^ W r l — M| + iM z T z 

where g is the charged current electro- weak coupling constant, Pl is the usual spin 
projector, 9 W is the Weinberg angle and r M represent the 4-momentum components of 
the Z boson, u and v are the spinors associated with va and vb, their functional form 
can be found in [19]. r 2 is at most the order of p 4 /M Pl and hence for any incoming 
momenta p < 10 19 ' 5 eV r 2 Mf . After some brute force algebra, we end with 

Mab x miL^ey iQEvEplEqEq ' x F{x) ' (10) 

where \ is the angle between incoming and outgoing va and F(x) is a complicated 
function of X- 

% We notice here that according to our eq. (jj) neutrino splitting is forbidden for superluminal neutrinos 
in the case n = 1 studied e.g. in [51]. Although the model presented in [53] does not have an EFT 
description, a simple threshold computation under the assumption of energy-momentum conservation 
and equal values of masses and LV parameters for neutrinos and antineutrinos shows that if neutrino 
splitting is possible also in that framework, values ~ — 10~ 51 can be probed by ~ 10 19 eV neutrinos. 
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The phase space, oI&ab is given by 

(27r) 3 2£ p , (27r) 3 2E q/ (2tt) 3 2£, 



We integrate over q' and we are left with 



Substituting Mab and d&AB back into (jSJ) we have 

J 5(E P -E' p -E q - E' q )d 3 p'd 3 qF 2 . (13) 



9 4 

AB 



16(2vr) 5 M|cos 4 ^ 

We now turn to estimating the size of the remaining integral well above threshold, 
as UHE neutrinos are far in excess of the neutrino splitting threshold around 20 TeV. 
Temporarily, we will set r\ v = 1 and re-insert it at the end. We first assume that the 
opening angle for all three neutrinos is roughly equal, as is true over most of phase space. 
In this case, F(x) can De approximated as (1 — cos x)- The opening angle x vanishes as 
Mpi oo, and is small even for UHE neutrinos. F(x) then reduces to F(x) = X 2 /2- 
Indeed, x is given by the typical transverse momenta p± divided by the longitudinal 
momenta p». The characteristic size of p± can be estimated by energy conservation, 
recognizing that when energy is conserved any "excess" energy from LV goes into the 
energy needed to create transverse momenta. Formally this implies that 

^ 7,' (14) 

i.e. p± ~ pn/Mpi- p|| itself for any particle can range from almost to almost the initial 
energy E p well above threshold. The available phase space volume for the remaining 
outgoing particles can therefore be approximated (although somewhat overestimated) 
as a cylinder with length E p and radius E 2 /Mpi, which has a total volume of nE^/M^. 
Similarly, x can De estimated as x ~ P\\/Mp\ ~ E p /Mp\. The 5-function in energy 
simply removes one factor of energy from our equation. Putting these expressions back 
into ffTBl we find 



(15) 



16(2tt) 5 M|cos 4 ^ E p \M 2 J 4M 4 , ' 
or, simplified 

r " ~ Af 3 ' (16) 

where we used G F j\f2 = g 2 /{%M 2 z cos 2 8 w ). If we now want to restore the r\ v factor, we 
notice that we must put one r\ v for each Mp 2 . In the end, the total rate is computed 
by adding the partial rates for each neutrino flavor, which is equivalent to multiply by 
3. Hence, we have 

r 3G 2 FV tE^ 



Possible cosmogenic neutrino constraints on Planck-scale Lorentz violation 



10 



This width can be turned into a decay length as 

L '" = ft ~ 17 X 10 " MPC ( lO^ev) U ' (18) 
This makes clear we need to push the required energies above 10 185 eV (with rj v = 1) 
for the rate to be appreciable. As a final remark, we notice that the decay length in 
Eq. fj!8j) strongly depends on both the energy and r\ v . Therefore, the error about its 
actual magnitude we might have made in our estimate will reflect in very small errors 
in the determination of the energy at which LV effects start to be relevant as well as of 
the constraint on r\ v . 

4-1.2. Z boson resonance At such high energies the Z could be real - i.e. there is a 
resonance in the matrix element. Even in this regime, however, the neutrino decay 
time can be computed easily, as the only hypothesis one has to relax is that the Z 4- 
momentum r satisfies r 2 <C M§. The magnitude of r 2 can be easily computed exploiting 
the kinematic equations. We obtain r 2 = 16/27i] l/ E A /Mp l . The final decay length is then 



L uup ~ 1.7 x 10- Mpc r,^ ( ] x 



27 Vv M^M 2 z J \M Z 



.(19) 



A comparison between the two different decay lengths, Eq. ( TT8]) and Eq. ( TT9|) . can be 
found in Fig. [JJ The Z resonance is hit at E A = 27/16 (rj^M^M^) . We notice that even 




log lo (E/eV) 

Figure 1. Comparison between computations of the decay length without 
(Eq. (fT5]) /red dashed line), and with (Eq. (TTTJ]) /black solid line) the Z boson resonance. 

though the two computations lead to very different results above the resonance, they 
will not lead to any appreciable effect in the neutrino spectra, as at such energies the 
decay lengths are anyway much smaller than the propagation distance of cosmological 
neutrinos. 
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4-2. A non-optimal case: flavor dependent LV 

In the flavor blind case above every mass eigenstate undergoes splitting. It is easy 
to construct a scenario where only one mass eigenstate decays - one eigenstate has a 
positive 7] LV coefficient while the other two have no LV. While this is unnatural in 
some sense, it serves as a nice example to illuminate a method to easily hide LV. 

The neutrino spectrum produced by the GZK process has a distinct distribution in 
flavor states of 1/3 v e and 2/3 v^. Therefore, in order to calculate a spectrum with LV 
we need to convert this flavor spectrum into a mass spectrum, i.e. we need to choose a 
neutrino mixing model. For our purposes we choose tribimaximal mixing [731 ES], which 
satisfies current experimental constraints and provides a simple mixing matrix. With 
tribimaximal mixing the GZK neutrino spectrum has equal distribution over all mass 
eigenstates. 

Now that we have spelled out out test models we shall proceed discussing the results 
for both scenarios. 

5. Results 

The net effect of neutrino splitting is very simple. It kills one neutrino of energy E 
and creates 3 neutrinos of average energy E/3, provided the parent neutrino is above 
threshold and has a reasonable life time. As we showed in sec. I4.1.1[ for a life time 
shorter than the age of the Universe, the neutrino energy has to be above 10 18 ' 5 eV, 
i.e. we need to probe UHE neutrinos. 

The effects of neutrino splitting on the UHE neutrino spectrum are twofold and 
can be understood qualitatively as follows. 

Flux suppression at UH energies The splitting is effectively an energy loss process 
for UHE neutrinos. If the rate is sufficiently high, the energy loss length can be 
below 1 Mpc. Let us call E(rj v ) the energy at which this happens. Then, being 
GZK neutrinos produced mainly at distances larger than 1 Mpc, we do not expect 
any neutrino to be detected at Earth with E > E. 

The mere observation of neutrinos up to a certain energy E ^ s would imply a 
constraint, according to Eq. IT51 



Flux enhancement at sub-UH energies Neutrinos lose energy by producing lower 
energy neutrinos. Eventually these neutrinos will become stable, either because 
their energy is below threshold, or because their lifetime is larger than their 
propagation time. Accordingly, we expect an enhancement of the neutrino flux 
at energies below few x 10 18 eV. 

In spite of being qualitatively straightforward, however, this analysis is not powerful 
enough to provide us with constraints on LV. In order to obtain meaningful constraints, 




(20) 
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we have to resort to full MonteCarlo simulations of the UHECR propagation from 
sources to the Earth. 

We simulated then the propagation of UHECR protons in the Inter Galactic 
Medium using the Monte Carlo package CRPropa [65], suitably modified to take into 
account LV in the neutrino sector. The simulation parameters are the following: we 
simulated unidimensional UHECR proton propagation, with source energy spectrum 
dN/dE oc E~ 2 ' 2 , from a spatially uniform distribution of sources located at redshift 
z < 3 according to the Waxman & Bahcall (WB) distribution used in [66]. The injection 
proton spectrum was tuned to fit AUGER data [59|. 

5.1. Results for flavor bind LV scenario 

Figure [2] shows the outcome of the simulations for different values of the LV parameter r\ v 
in the best case scenario, together with experimental sensitivities from some existing and 
planned observatories, as well as the Waxman & Bahcall bound [68l 169] for reference^ 
Results are in agreement with qualitative expectations previously discussed: Above a 
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Figure 2. Evolution of the predicted LV neutrino spectra varying f] v in the "best case 
scenario" . Sensitivities of main UHE neutrino operating and planned experiments are 
shown, as found in [2U LHJ EZ] • The Waxman & Bahcall limit 68, 69] in the interesting 
energy range is shown for reference. 

+ This limit is in fact an estimate of the neutrino luminosity of sources of UHE Cosmic Rays and 
7-rays, in the hypothesis that the sources are optically thin to the escape of UHE particles and that 
both 7-rays and neutrinos are originated from UHECR interactions with radiation backgrounds. It is 
worth mentioning that this bound might be strongly affected by QG effects, as shown in |70) . 
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certain energy the neutrino spectrum displays a sharp cut off, while at lower energy a 
peak appears. The peak position and strength, as well as the position of the cut off, 
depend on the value of the LV parameter r\ v . If 77^ = 1 the peak is found at < 10 18 eV 
and overwhelms the LI spectrum by a factor of roughly 10, while the cut off energy is 
at ~ 10 18 ' 7 eV. Also the evolution of the neutrino spectra with r\ v can be seen in the 
same Fig. |2j 

Although at present it is not possible to draw firm conclusions for the constraints 
on LV in the neutrino sector, we notice that future experiments, if not already AUGER, 
will be able to probe the fluxes we predict, hence they will be able to cast constraints 
on r\ v . In the future, the ARIANNA experiment [50J will be able to probe fluxes down 
to ~ 3 eVcm~ 2 s^ 1 sr^ 1 in the energy range 10 17 -r- 10 20 eV in the full configuration 
[5T] . If such a sensitivity will be achieved, a constraint of order r\ v < 10~ 4 will be cast, 
according to Fig. [2j Moreover, the IceRay experiment [52] planned at the South pole 
is expected to observe roughly 4 neutrino events per year of data taking for the WB 
model preferentially studied here, in the range 10 17 10 19 - 5 eV [52J. Hence, constraints 
of order rj < 10 -3 are expected after few years of data taking by this experiment. We 
finally notice here that by exploiting this strategy nothing can be said about the case 
Vu < o. 

5.2. Results for flavor dependent LV scenario 

If only one mass eigenstate decays, the energy momentum conservation equations, rates, 
and lifetimes all remain the same. We therefore can apply our MonteCarlo spectrum 
for this scenario to just any one of the mass eigenstates (given our assumption of 
tribimaximal mixing no neutrino mass eigenstate is preferred). The only difference 
is that only 1/3 of the total flux undergoes neutrino splitting. The resulting spectrum 
is shown in Fig. [3] Note that even for rj — 1 the departure from the LI spectrum is 
significantly reduced, as one would expect, with the magnitude of the deviation at the 
peak being lowered to a ~ 50% excess over the LI flux. Hence, with present sensitivity, 
absence of a signal in the neutrino spectrum does not mean that LV is ruled out for 
neutrinos, even for 0(1) positive coefficients. We notice that the cutoff feature present 
in the flavor blind case has now effectively disappeared. Then, observations of a UHE 
neutrino flux up to some maximal energy do not allow to draw conclusions about LV in 
this case scenario. 

We now return to the question of negative coefficients and neutrino splitting. So 
far, in both our flavor blind and flavor dependent scenarios all the coefficients have been 
positive. As argued in section HI neutrinos can split even with negative coefficients when 
the relation n UB < r] UA — K(p) holds. However, one can immediately see that therefore 
at least one of the mass eigenstates is always stable if the coefficients are negative. This 
implies that no equivalent of the "best case scenario" (all the mass eigenstates decay) is 
possible for negative rj u i and that, consequently, the deviation from the LI spectrum will 
be suppressed by at least 1/3. Therefore cosmogenic neutrinos are not effective probes 
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Figure 3. Observed UHE neutrino spectra for different values of the LV coefficient 
of the only mass eigenstate undergoing splitting in the non-optimal scenario. The flux 
suppression at UH energies has effectively disappeared (for what concerns observational 
relevance) and the excess at lower energies is significant only for O(l) LV coefficient. 



of the 7] uI < region of parameter space. 
5.3. Dependence on model uncertainties 

Predictions of cosmogenic neutrino spectra are known to be plagued by several 
uncertainties, as they depend strongly on the evolution of UHECR sources with redshift, 
which cannot be directly probed with UHECRs due to the GZK attenuation effect. As 
a consequence, also our predictions for the UHE neutrino flux in the presence of LV can 
be affected. 

We check this by simulating different UHECR sources evolution models for the "best 
case scenario" . In Fig. H] we show the results for the model outlined in [SB] , which is 
expected to give the largest flux. As a result, the presence and the strength of the bump 
feature are demonstrated to be weakly dependent on the underlying UHECR source 
distribution, with the choice of the source distribution affecting by a factor of 2 the 
height of the bump. This was somewhat expected from our qualitative considerations: 
as long as UHECR sources are located at distances much larger than 1 Mpc from 
Earth, the magnitude of the bump depends only on how many neutrinos were originally 
produced at energy larger than the LV cutoff energy. On the other hand, the cutoff 
feature is definitely model independent, as long as the unknown maximal energy at 
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Figure 4. Evolution of the predicted LV neutrino spectra varying r/^ for the maximal 
model 1 56 . 



which UHECRs are accelerated is enough to produce neutrinos at energy larger than 



6 x 10 18 ■ ?T 4/13 eV. 



6. Other processes 

As we mentioned in the previous sections, there are other processes, besides splitting, 
leading to neutrinos losing energy in the inter-galactic medium. In particular, the 
hadronic decay modes of the Z° could lead to the production of UHE protons, thereby 
mimicking a "Z-burst" effect [7T] . Although this model was conceived for other purposes 
than constraining LV, and has now lost most of its attraction due to experimental results 
on proton spectra above the GZK threshold, it is notable that the same mechanism could 
in principle help constrain LV in the neutrino sector. However, we can argue with a 
very simple argument that this is not the case. 

Let us consider the LI neutrino flux, which corresponds in Fig. [2] to the rj = case. 
The minimal process leading to proton production in this context is v — > vZ* — > upp, 
whose threshold energy, if t] u = 1, is at ~ 10 18 ' 7 eV. Processes involving more particles 
are naturally suppressed, however they represent the majority of the allowed channels 
and hence they might be relevant to our case. Even if we do not know the energy 
spectrum of the produced protons, we know that their mean energy is roughly 1/80 times 
the energy of the parent neutrino [7T] . We can then compute the maximal expected flux 
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of UHECRs produced by this Z-burst-like mechanism by converting all the neutrinos 
with energy E v above threshold into 2 protons of energy 1/80 x E v . The result is 
shown in Fig. |5j This mechanism would indeed contribute to the UHECR spectrum at 




log io (E/eV) 

Figure 5. Observed UHE proton spectrum compared with predictions for just primary 
(yellow) and "Z-burst"-like originated protons (blue). The latter might produce an 
excess in the UHECR spectrum at 10 18 eV, no larger than 5%, and hence below the 
systematic uncertainties of UHECR experiments. Note that, in order to speed up 
computation, we switched off pair production, which is known to play an important 
role in reproducing the ankle feature 72 . This is why our UHECR spectrum is above 
data between 10 185 eV and 10 19 eV. 

10 18 eV, producing an excess no larger than 5%, which however is below the systematic 
uncertainties of UHECR experiments. Therefore, no constraints on LV can be implied 
using this technique. 

Electron/positron pair production would also contribute to the cutoff in the 
neutrino spectrum. The threshold for pair production of electrons is at roughly 
V^vMpi = 10 17 eV. At energies near the existing neutrino splitting cutoff, 10 18,5 eV we 
are well above threshold. The mass of the electrons is therefore irrelevant in calculating 
the decay rate (just as the neutrino mass was) and hence the electron pair production 
process has nearly the same rate as the neutrino splitting process. If neutrinos produce 
electron/positron pairs (which after all depends on LV in the electron sector) it can only 
increase the steepness of the cutoff. It will not contribute to the excess neutrino flux 
at 10 17 ' 5 eV and hence one needs to probe energies one order of magnitude larger than 
10 17 ' 5 eV to derive constraints exploiting the presence or the absence of a cutoff feature. 
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7. Conclusions 

In this work we have investigated possible signals of higher dimension LV CPT even 
operators in the UHE neutrino spectrum, in particular the effect of the neutrino splitting 
on the UHE neutrino spectrum. This process provides a clean test as it does not involve 
other LV operators apart from the neutrinos' ones. In addition, since the dominant 
neutrino field is left-handed, there are no complications due to different LV for different 
chirality of fermion, as there are in the UHECR case, and one ends up with one LV 
coefficient for each neutrino mass eigenstate. In the flavor blind scenario, where every 
mass eigenstate undergoes splitting approximately at the same energy, there is both a 
precocious fall off of the neutrino flux at UHE as well as a significant excess in the UHE 
neutrino flux at energies as low as 10 175 eV. Noticeably, this kind of energies are well 
within reach of current and future UHECR experiments and about order of magnitude 
below those so far used for LV tests with hadronic and electromagnetic UHECRs. 

According to our study, existing or planned UHE neutrino experiments have the 
potential to probe LV in the neutrino sector for coefficients r] > 10 -4 . However, we have 
discovered a serious difficulty with deriving constraints in the absence of a positive LV 
signal event. The distribution in mass eigenstates is roughly equal for UHE neutrinos 
in realistic mixing scenarios. Although it might seem somewhat unnatural that LV has 
different effects on different mass states of the same particle field, if this is the case then 
it is possible that LV can exist/be strong for one mass eigenstate yet be almost invisible 
for UHE neutrino detectors. In particular, the observation of a neutrino flux up to some 
maximal energy does not imply a firm conclusion on LV, at least with present accuracy. 
However, let us note that if a neutrino is ever detected at energy E u , we can rule out 
a flavor blind r\ v > (E^/IO 18 - 8 eV)~ 13 / 4 according to Eq. (120]) . On the other hand, the 
bump feature can lead to constraints on r] u > 0, as this bump should be observable in 
any such scenario of LV in the neutrino sector. In this case, to obtain a 0(1) constraint 
on LV requires at least a 50% accuracy in the determination of the neutrino spectrum 
in the energy range 10 17 -j- 10 19 eV, which can be achieved by future experiments. 
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